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We use all-order relativistic many-body perturbation theory to study 5s^nl configurations of In I 
and Sn II. Energies, El-amplitudes, and hyperfine constants are calculated using all-order method, 
which accounts for single and double excitations of the Dirac-Fock wave functions. A comprehensive 
review of experimental and theoretical studies of In I and Sn II properties is given. Our results are 
compared with other studies were available. 
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I. INTRODUCTION 

In this work, we present a systematic calculation of 
various In I and Sn II atomic properties and study the 
importance of the high-order correlation corrections to 
those properties using relativistic all-order method. Pre- 
viously these atoms have been studied in a number of 
experimental and theoretical papers. First theoretical 
studies were published 30 years ago by Migdalek 
They used relativistic semiempirical method including 
exchange to calculate the oscillator strengths in In I for 
the 5s^5pj - hs^Qpj - 5s^7isi/2, 5s^5pj - 5s^ndj, 

and 5s^6si/2 - 5s^npj transitions. 

Later, the oscillator strengths determined from single- 
configuration relativistic Hartree-Fock (RHF) calcula- 
tions were reported by Migdalek and Baylis \2\ for the 
lowest 5s^5pj-5s^6si/2 and bs^5pj-5s^5dji transitions. 
A quantum defect theory was used by Gruzdev and 
Afanaseva [3] to calculate oscillator strengths / averaged 
over j in neutral indium. Configuration interaction gf 
values for transitions between the 5s^6si/2, 5s^ndj, and 
bs^npj (with 71 = 5, 6) states were reported for the in- 
dium isoelectronic sequence up to Ba VIII in Ref. [1]. 
A self-consistent-field method was used to generate one- 
electron orbitals. The method used in 0] included rela- 
tivistic effects albeit in an approximate way, and the con- 
figuration interaction scheme accounts for correlation ef- 
fects [^] . Hartree-Fock calculations including relativistic 
corrections and configuration interaction in an intermedi- 
ate coupling scheme were carried out in Ref. [s] to analyze 
the spectrum of Sn II. Transition probabilities for 36 lines 
of Sn II arising from the 5s^ns, bs^np, bs^nd, bs^nf , and 
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5s5p^ configurations of Sn II were evaluated in [s^] using 
the Cowan code. Radiative transition probabilities and 
oscillator strengths for 164 lines arising from the bs^ns, 
bs^np, bs^nd, bs^nf, bs^ng, and bsbp^ configurations of 
Sn II were calculated recently by Alonso-Medina et al. in 
[9] . These values were obtained in intermediate coupling 
(IC) using ab initio relativistic Hartree-Fock (HFR) cal- 
culations. The standard method of least square fitting of 
experimental energy levels by means of computer codes 
from Cowan was used @ to calculate IC transition rates. 
Recently, energies of the bs^bpj, bsbp^, 5s^6si/2, bs'^bdj, 
and bs^Gpj states in Sn II were evaluated by Dzuba and 
Flambaum in [t^ using many-body perturbation theory 
(MBPT). It was underlined that correlations and rela- 
tivistic corrections were important. The screening of the 
Coulomb interaction and hole-particle interaction was in- 
cluded in all orders of the MBPT [7]. 

The experimental study of atomic lifetimes in gal- 
lium, indium, and thallium was carried by Andersen and 
S0rensen Q using beam-foil technique. Results for the 
bs^6s and bs^nd (n = 5-7) levels in In I were given 
in Q. Lifetimes of the bs^ns and bs^nd {n < = 20) 
states in indium measured using pulsed laser excitation 
of an atomic beam were reported by Jonsson et al. in 
0. Determination of radiative lifetimes of the bs^6p, 
bs'^ns, and bs^nd {n < = 10) levels in In I using a pulsed 
laser was presented in Refs. |10l.[T]|. The atoms were ex- 
cited in an atomic beam, with a nitrogen-laser-pumped 
dye laser. The fiuorescence decay from the atoms was 
observed by a fast photomultiplier [l^, [IH. The op- 
tical emission from a laser produced plasma generated 
by 1064 nm irradiation of Sn/Pb alloy targets at a flux 
of 2x10^*^ Wcm^^ was recorded and analyzed between 
200 and 700 nm [5]. Experimental transition probabili- 
ties for 36 lines of Sn II arising from the bs^ns, bs^np, 
bs^nd, bs^nf, and bsbp^ configurations of Sn II were de- 
termined by Alonso-Medina et al. in Q . Lifetime mea- 
surements for levels arising from the bs^bd and bs^Af 
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configurations in Sn II were presented by Schectman et 
al. in ^121] . These measurements utilized the Univer- 
sity of Toledo Heavy Ion Accelerator Beam Foil Facility. 
The results were discussed in the context of interpret- 
ing vacuum ultraviolet absorption spectra observed with 
the Goddard High Resolution Spectrograph on board the 
Hubble Space Telescope 

Hyperfine structure of the hs^bpj states of In^^^ and 
In-'^^'^ were measured by the magnetic-resonance method 
An atomic beam irradiated by a narrow band 
dye laser was used in [l^ to observe resonance fluo- 
rescence in free indium atom. From the resonance fre- 
quencies, values for the hyperfine structure of the 5s^6s 
level in In^^^ and In^^^ were derived by Neijzen and 
Donszelmann in [l^. The spin- forbidden 5s^5p 
- 5s5p^ In^^^ transition was analyzed and absolute 
wavelengths, hyperfine constants A and B, as well as im- 
proved energy level values were reported by Karlsson and 
Litzen in ^16] . 

A high-resolution study of the A = 451.1 nm transi- 
tion in In I using CW dye laser was reported by Zaal 
et al. in [l^]. The blue dye laser set-up was tested on 
the 5s^5p3/2-5s^6si /2 A — 451.1 nm transition in natural 



indium [17[ . Proposal for laser on self-terminating tran- 
sition in blue spectral range on indium atom transition at 
451.1 nm was presented recently by Riyves et al. in [l8| . 
The spectroscopy of dense In vapor was studied recently 
via resonant pulsed laser excitation at A = 410.13 nm 
(the 5s^5pi/2-5s^6si/2 transition) (loj . 

In this paper, we conduct both relativistic many- 
body perturbation theory (RMBPT) and all-order single- 
double (SD) calculations of In I and Sn II properties. 
Such calculations permit one to investigate convergence 
of perturbation theory and estimate the uncertainty of 
theoretical predictions. We evaluate reduced matrix ele- 
ments, oscillator strengths, and transition rates for pos- 
sible 5s^n/ — 5s^n'l' electric-dipole transitions in In I and 
Sn II and calculate the lifetimes of the corresponding 
levels. Our results are__compared with theoretical results 
from Refs. [1 
Refs. 



I , Q and with measurements from 
3^ Inf in In I and Refs. [1 HI in Sn II. 
We also calculate hyperfine constants A for the b-s^npj 



in- 



5s^nsi/2 [n = 6 — 9), and hs^ndj (n = 5- 



states in ^^^In using the relativistic MBPT and SD all- 
order methods. Where possible, we compare our results 
with the measurements from Refs. [l^ [l4. ITsj. 

We consider the three-electron system 
[Ni]4s^4p^4d^°5s^nZ in In I and Sn II as a one-electron 
nl system with [Ni]4s^4p^4(i^°5s^ core. Recently, the 
relativistic all-order method was used to evaluate the 
excitation energies, oscillator strengths, transition rates, 
and lifetimes in Ga I [1^ as well as in Tl I and Tl-like 
Pb dij. The [Ni]4s^nZ states in Ga I were treated in 
[20| as the nl one-electron system with [Ni]4s^ core 
and [Xe]Af^^bd'^%s'^nl states in Tl I and Pb II were 
evaluated in 21:] as the nl one-electron system with 
[Xe]4/i45di"6s2 core. 

To summarize, this work presents both a systematic 



calculation of various properties of In I and Sn II, and a 
study of the importance of the high-order correlation cor- 
rections to these properties. We conclude that all-order 
SD method, in general, produce more accurate values 
than the third order MBPT and can be used for the accu- 
rate calculation of In and Sn"*" properties. By comparing 
the all-order and third-order MBPT results, we were able 
to study the relative importance of the correlation cor- 
rections for different properties and single out the cases 
where the treatment of In as a three-particle system may 
be important, i.e. the cases where significant discrepan- 
cies between theory and experiment persist even for the 
all-order calculations. The development of the all-order 
approach that is capable to fully treat In or Sn+ as a 
three-particle system is a difficult problem [2^, H^, [l^l 
and the initial studies of the applicability of the all-order 
method to such systems may be useful. We find that 
the all-order SD method works relatively well for In I 
even without explicit consideration of the three-particle 
states. For Sn II, the convergence of MBPT expansion 
is worse than for In, particularly for d-wave, where SD 
equations diverge. That is caused by the strong inter- 
action between 5s^nd configurations and low-lying 5s5p^ 
configuration, which corresponds to the excitation from 
the core. 

In the next section, we briefly review the RMBPT the- 
ory and all-order SD method for the calculation of atomic 
properties of the atoms with one unpaired electron. The 
energies are given in Table IH Extension of the theory 
to one-electron matrix elements is discussed in Sec. HI. 
Our results for El transitions are listed in Tables UIl - fVl 
Calculated and experimental lifetimes for In are given in 
Table Magnetic hyperfine structure of In is discussed 
in Sec. IV and results are summarized in Table IVTIl 



II. ENERGIES OF In I AND Sn II 

We start from the "no-pair" Hamiltonian [H 
second quantization form 

H = Ho + Vi, 

Hq = ^ Sia'lai , 
i 

ijkl 



m 



the 

(1) 
(2) 

(3) 



where negative energy (positron) states are excluded 
from the sums; Si are eigenvalues of the one-electron DF 
equations with a frozen core, and gijki is the Coulomb 
two-particle matrix element. 

Considering neutral In as a one-electron system we use 
y^-i DF potential [Ni]4s24p^4rfi°5s2 to calculate DF 
orbitals and energies Si. There are a number of advan- 
tages associated with this potential, includin g a greatly 
reduced number of the Goldstone diagrams [27|, which 
leads to important simplifications in calculation. For ex- 
ample, when considering the total energy of different va- 
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TABLE I: Valence energies in different approximations for In I and Sn II in cm ^. We calculate zeroth-order (DF), single- 
double Coulomb correction E^^ , and the part of third-order -Eg^tra which is not included in the E^^ . Breit corrections 
B^"' are calculated in first and second orders. The sum of these five terms Ef^ is compared with experimental energies 
EmsT SE^^ = Ef^, - EmsT- The differences JS^^) and 5E^^^ between total energies = E^°> + E^^'> + B^'^^ + B^^\ 

-E'tot = -B tot + -E*^^) and experimental energies JSnist .251 are given for comparison. 



nlj 


^DF 


^SD 


^(3) 
extra 


5(1) 


5(2) 


tot 


-Enist 




5E^'^ 














In I 












5pi/2 


-41507 


-5554 


913 


105 


-146 


-46189 


-46670 


-2163 


547 


481 


5p3/2 


-39506 


-5378 


912 


73 


-132 


-44031 


-44457 


-2040 


557 


426 


5^3/2 


-12390 


-1350 


161 


2 


-4 


-13581 


-13778 


256 


493 


197 


5^5/2 


-12374 


-1337 


160 


1 


-4 


-13554 


-13755 


260 


493 


201 


6Sl/2 


-20572 


-2096 


232 


12 


-18 


-22442 


-22297 


-237 


249 


-145 


6pi/2 


-13979 


-964 


113 


13 


-16 


-14833 


-14853 


-150 


122 


20 


6P3/2 


-13719 


-919 


112 


10 


-16 


-14532 


-14555 


-155 


115 


23 


6^3/2 


-6955 


-558 


69 


1 


-2 


-7445 


-7809 


375 


479 


364 


64/2 


-6946 


-554 


68 


1 


-2 


-7433 


-7697 


277 


379 


264 


7Sl/2 


-9867 


-584 


72 


4 


-6 


-10381 


-10368 


-95 


64 


-13 


7pi/2 


-7488 


-349 


41 


5 


-6 


-7797 


-7809 


-56 


45 


12 


7p3/2 


-7388 


-335 


41 


4 


-6 


-7684 


-7697 


-58 


42 


13 


8Sl/2 


-5816 


-251 


32 


2 


-3 


-6036 


-6033 


-45 


27 


-3 


4/5/2 


-6863 


-118 


14 








-6967 


-6963 


-4 


9 


-4 


4/7/2 


-6863 


-118 


14 








-6967 


-6962 


-5 


8 


-5 


5/5/2 


-4393 


-67 


8 








-4452 


-4450 


-3 


5 


-2 


5/7/2 


-4393 


-67 


8 








-4452 


-4450 


-3 


5 


-2 


7^3/2 


-4441 


-285 


35 





-1 


-4692 


-4834 


148 


201 


142 


7^5/2 


-4436 


-283 


35 





-1 


-4685 


-4808 


129 


181 


123 


8pi/2 


-4687 


-168 


20 


2 


-3 


-4836 


-4843 


-28 


21 


7 


8P3/2 


-4638 


-162 


20 


2 


-3 


-4781 


-4789 


-28 


21 


8 


8(^3/2 


-3078 


-171 


21 





-1 


-3229 


-3334 


112 


144 


105 


8rf5/2 


-3075 


-170 


21 





-1 


-3225 


-3315 


98 


129 


90 


9Sl/2 


-3837 


-132 


17 


1 





-3951 


-3951 


-24 


15 















Sn II 












5pi/2 


-111452 


-6848 


1014 


206 


-233 


-117313 


-118017 


-2578 


712 


704 


5p3/2 


-107358 


-6719 


1032 


146 


-216 


-113115 


-113766 


-2459 


736 


651 


6Si/2 


-57995 


-3597 


467 


35 


-44 


-61134 


-61131 


-641 


435 


-3 


6pi/2 


-44483 


-2244 


271 


39 


-38 


-46455 


-46523 


-316 


314 


68 


6p3/2 


-43691 


-2133 


258 


28 


-37 


-45575 


-45640 


-297 


300 


65 


7Sl/2 


-30735 


-1230 


172 


14 


-17 


-31796 


-31737 


-343 


55 


-59 


4/5/2 


-27689 


-1189 


147 





-1 


-28732 


-28731 


-24 


206 


-1 


4/7/2 


-27691 


-1193 


147 





-1 


-28738 


-28725 


-34 


196 


-13 


7pi/2 


-25253 


-930 


115 


16 


-16 


-26068 


-26114 


-107 


153 


46 


7p3/2 


-24917 


-895 


109 


12 


-15 


-25706 


-25751 


-94 


152 


45 


8Sl/2 


-19133 


-629 


84 


7 


-9 


-19680 


-19615 


-162 


26 


-65 


5/5/2 


-17759 


-670 


85 





-1 


-18345 


-18358 


-13 


131 


13 


5/7/2 


-17761 


-674 


86 





-1 


-18350 


-18352 


-22 


121 


2 


8pi/2 


-16354 


-487 


60 


8 


-8 


-16781 


-16821 


-31 


102 


40 


8p3/2 


-16179 


-472 


57 


6 


-8 


-16596 


-16630 


-29 


97 


34 


9Sl/2 


-13070 


-301 


47 


4 





-13320 


-13337 


-88 


13 


17 
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lence states of a one-electron atom, that energy can be 
written as 



(4) 



where i?corc is the same for aU valence states v. The first- 
order correlation correction to valence removal energies 
vanishes for a V^~^ DF potential and the first nonvan- 
ishing corrections appear in the second order (28[ : 



-EE 

mn a 

+EE 



n ah 



Qavmn {^Qmnav ffmnva ) 
9nvha{,9ahnv Qahvn) 



(5) 



We use indexes a and b to label core states and m and n to 
designate any excited states. The second-order Coulonib- 

(2) (2) 

Breit contribution Bi, ' is obtained from the Ei, ' expres- 
sion (O by changing gijki gijki+bijki and keeping only 
terms that arc linear in bijki that is a two-particle matrix 
element of the Breit interaction [29] : 

B = — — [aia2 - ^ [aia2 - (Q;ifi2) (q;2?12)]] , (6) 

where ai is the Dirac matrix, ri2 — ri2/ri2, and a is the 
fine structure constant. The first-order Breit correction 

is B^^' = J2a [bvava - Kaav] ^ ~ J2a ^vaav, whcrC dirCCt 

term vanishes after summing over closed shells. 

Even though the number of Goldstone diagrams for 
the DF potential is much smaller than in general 

case, the third-order expression for energy correction still 
includes 52 terms. The corresponding formula for E'i^' 
was presented by Blundell et al. in Ref. (s^l, where 52 
terms were combined into 12 groups with distinct energy 
denominators: 



= E 



(3) 



E 



(3) 



(7) 



Expression ([7]) includes terms with one-, two-, three-, and 
four-particle sums over virtual states in addition to sums 
over core states. 

The all-order SD method was discussed previously in 
Refs. [2i|,[3l|,[33,[3l[3l,[3i,|3|. Briefly, we represent the 
wave function of the atom with one valence electron 
as = -^f^: 



SD 



E + 



E 

mnab 



(8) 



1 + + 



where is the lowest-order atomic wave function, which 
is taken to be the frozen-core DF wave function of a 
state V. The coupled equations for the single- {pm^ and 
Pma) and double-excitation coefficients Pmnva and Pmnab 
are obtained by substituting the wave function ^'^^ into 



the many-body Schrodinger equation, with Hamiltonian 
given by Eqs. ((T]-[3]). Note that we again start from y"^^ 
DF potential. The coupled equations for the excitation 
coefficients are solved iteratively. In the following sec- 
tions, the resulting excitation coefficients are used to 
evaluate hyperfine constants and transition matrix ele- 
ments. 

The valence energy E^^ is given by 



' ^ ^ UvamnPmnva ■ 
mna 



— ^ ^ UvavmPma^ ^ ^ IJabvmPmvab 
ma mab 

(9) 

This expression does not include a certain part of the 
third-order MBPT contribution. This part of the third- 
order contribution, £^^^^^^^.3^, is given in Ref. ^34'| and has 
to be calculated separately. We use our third-order en- 
ergy code to separate out Ey^l^^^.^ and add it to the E^^ . 

We drop the index v in the eP, Ei^\ and E^^ desi gna- 
tions in the text and tables below. 

We use B-splines 37] to generate a basis set of DF 
wave functions for the calculations of MBPT and all- 
order expressions. Typically, we use 40 or 50 splines of 
order fc = 7 or 9, respectively, for each partial wave (see 
below for more details). Basis orbitals for In I and Sn II 
are constrained to cavities of radii i? = 95 a.u. and R — 
85 a.u. respectively. The cavity radii are chosen large 
enough to accommodate all orbitals considered in this 
paper and small enough for 50 splines to approximate 
inner-shell DF wave functions with good precision. 

Results of our all-order SD calculations of energies for 
the lowest states of neutral In and In-like Sn ion are given 
in Table |T1 Our final answer Ef,^ also includes the part 
of the third-order energies omitted in the SD calculation 
as well as the first-order Breit correction B'^^^ and 
the second-order Coulomb-Breit i?^^) correction. Theo- 
retical values are compared with the recommended values 
EmST from the National Institute of Standards and Tech- 
nology database [111, SE^^ = Ef^ - i^NiST- For com- 
parison, we also give the differences between the second- 
order and third-order MBPT calculations and experi- 
mental values in columns labelled (5£''^' and 6E'-^\ In 
Sn II the all-order SD equations for d-wave do not con- 
verge and we exclude d-orbitals of Sn II from Table HI 

The largest correlation contribution to the valence en- 
ergy comes from the second-order term, i?*^^-* . As we have 
discussed above, this term is simple to calculate in com- 
parison with E^^") and E^^ terms. Thus, we calculate 
E'--^^ with better accuracy than E^'^'' and E^^ . To in- 
crease the accuracy of the i?^^-* calculations, we use 50 
splines of order k = 9 for each partial wave and include 
partial waves up to Zmax = 10. Then, the final value is ex- 
trapolated to account for contributions from higher par- 
tial waves (see, for example, Refs. [38, 39j).We estimate 
the numerical uncertainty of E^-^^^ caused by incomplete- 
ness of the basis set to be approximately 10 cm""'^ or less, 
depending on the valence state. 

Owing to the numerical complexity of the E^^ calcula- 
tion, we use Imax = 6 and 40 splines of order k = 7. As we 
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TABLE II: Wavelengths A (A), transition rates Aj. (s oscillator strengths (/), and line strengths S (a.u.) for transitions in 
In I calculated in all-order perturbation theory. Numbers in brackets represent powers of 10. 



Transition 


A 




J 










Transition 


A 




J 








5pi/2 


6Sl/2 


4153 


5.15[7] 


1.33 


-1] 


3.64 


0] 


4/5/2 


84/2 


26532 


3.32[5] 


3.51 


-2] 


1.84[ 


1 1 


5pi/2 


5(^3/2 


3045 


1.30 [8] 


3.61 


-1] 


7.24 


0] 


4/7/2 


74/2 


42937 


1.65 [5] 


3.42 


-2] 


3.86[ 


1 1 


5pi/2 


7Sl/2 


2792 


1.37[7] 


1.60 


-2] 


2.93 


-1] 


4/7/2 


84/2 


26532 


6.65[4] 


5.26 


-3] 


3.68[ 


OJ 


5pi/2 


6^3/2 


2578 


3.o3[7J 


7.24 


-2] 


1.23 


0] 


8pi/2 


74/2 


529101 


4.67[3] 


3.92 


-1] 


1.37[ 


3J 


5Pl/2 


8Sl/2 


2493 


5.81 [6] 


5.41 


-3] 


8.88 


-2] 


8Pl/2 


84/2 


61275 


3. 79 [5] 


4.27 


-1] 


1.72| 


2J 


5pi/2 


7d3/2 


2410 


1.47[7] 


2.56 


-2] 


4.06 


-1] 


8P3/2 


74/2 


740741 


3. 40 [2] 


2.80 


-2] 


2.73| 


2J 


5pi/2 


8^3/2 


2329 


7. 49 [6] 


1.22 


-2] 


1.87 


-1] 


8p3/2 


74/2 


704225 


2.38[3J 


2.66 


-1] 


2.47[ 


3J 


5p3/2 


6Si/2 


4576 


9.05[7] 


1.42 


-1] 


8.56 


0] 


8p3/2 


84/2 


63371 


r /1 1 
8.38[4J 


5.05 


-2] 


4.21[ 


1 1 

IJ 


5p3/2 


5(^3/2 


3266 


2.49[7J 


3.98 


-2] 


1.71 


0] 


8p3/2 


84/2 


63211 


4.95 [5] 


4.44 


-1] 


3.70[ 


Ol 

2J 


5p3/2 


5(^5/2 


3264 


1.47[8J 


3.53 


-1] 


1.52 


1] 


5,f5/2 


84/2 


79681 


1.22[5] 


7.75 


-2] 


1.22[ 


Ol 

2J 


5P3/2 


7Sl/2 


2977 


2.31[7] 


1.53 


-2] 


6.02 


-1] 


5,f5/2 


84/2 


79428 


5. 70 [3] 


5.39 


-3] 


8.46[ 


nl 

OJ 


5P3/2 


6^3/2 


2734 


6.88[6J 


7.71 


-3] 


2.78 


-1] 


5/7/2 


84/2 


79428 


1.14[5] 


8.09 


-2] 


1.69| 


ol 

2J 


53?3/2 


6^5/2 


2734 


4.05[7] 


6.81 


-2] 


2.45 


0] 


6Sl/2 


6pi/2 


13669 


1.43[7] 


4.02 


-1] 


3.61| 


-1 1 

IJ 


5P3/2 


8Si/2 


2639 


9.72[6J 


5.07 


-3] 


1.76 


-1] 


6Si/2 


6P3/2 


13146 


1.57[7] 


8.13 


-1] 


7.03[ 


1 1 

IJ 


5p3/2 


7rf3/2 


2547 


2.79[6J 


2.71 


-3] 


9.08 


-2] 


6Si/2 


7pi/2 


7002 


1.40 [6] 


1.03 


-2] 


4.75[ 


1 1 

-IJ 


53?3/2 


7d5/2 


2546 


1.63 [7] 


2.38 


-2] 


7.97 


-1] 


6Si/2 


7p3/2 


6949 


1.96 [6] 


2.84 


-2] 


1.30[ 


OJ 


5P3/2 


8^3/2 


2456 


1.43 [6] 


1.29 


-3] 


4.18 


-2] 


6si/2 


8pi/2 


5806 


4.07[5] 


2.06 


-3] 


7.86[ 


ol 

-2J 


5P3/2 


8^5/2 


2456 


8.32[6J 


1.13 


-2] 


3.65 


-1] 


6si/2 


8P3/2 


5787 


6. 40 [5] 


6.42 


-3] 


2.45[ 


1 1 

-IJ 


6pi/2 


5d3/2 


69156 


1.57[5] 


2.25 


-1] 


1.03 


2] 


54/2 


7Pl/2 


18116 


7.58[5] 


1.86 


-2] 


4.45| 


nl 

OJ 


6pi/2 


7Si/2 


22594 


3. 48 [6] 


2.66 


-1] 


3.96 


1] 


54/2 


7p3/2 


17765 


6.13[4J 


2.90 


-3] 


6.78| 


1 1 

-ij 


6pi/2 


6^3/2 


13512 


7.35[6] 


4.02 


-1] 


3.58 


1] 


54/2 


4/5/2 


15798 


1.32 [7] 


7.43 


-1] 


1.55[ 


ol 

2J 


6pi/2 


8Sl/2 


11463 


1.13[6] 


2.22 


-2] 


1.68 


0] 


54/2 


8pi/2 


11816 


3.07[5] 


3.21 


-3] 


5.00| 


1 1 

-IJ 


6P1/2 


7d3/2 


9903 


3.79[6J 


1.12 


-1] 


7.27 


0] 


54/2 


8p3/2 


11741 


2. 46 [4] 


5.09 


-4] 


7.87| 


ol 

-2J 


6pi/2 


8(^3/2 


8665 


2.19[6] 


4.92 


-2] 


2.81 


0] 


54/2 


5/5/2 


11312 


5. 54 [6] 


1.59 


-1] 


2.37[ 


1 1 

IJ 


6P3/2 


5^3/2 


86580 


1.60 [4] 


1.80 


-2] 


2.05 


1] 


54/2 


7p3/2 


17838 


5.67[5] 


1.80 


-2] 


6.36[ 


nl 

OJ 


6P3/2 


5(^6/2 


84890 


1.03[5] 


1.66 


-1] 


1.86 


2] 


54/2 


4/5/2 


15855 


9.46[5] 


3.56 


-2] 


1.12[ 


1 1 

ij 


6P3/2 


7Sl/2 


24184 


6.3o[6J 


2.79 


-1] 


8.88 


1] 


54/2 


4/7/2 


15855 


1.42[7J 


7.13 


-1] 


2.23| 


ol 

2J 


6p3/2 


6^3/2 


14065 


1.54 [6] 


4.56 


-2] 


8.45 


0] 


54/2 


8p3/2 


11773 


2.29[5] 


3.17 


-3] 


7.36[ 


1 1 

-IJ 


6p3/2 


6^5/2 


i4U4o 


9.iz[Dj 


4.U4 


1 1 
-ij 


( .45 


-1 1 
ij 


54/2 


5/5/2 




3.93[oJ 


'yen 


-3J 


1. ( U| 


nl 

OJ 


6P3/2 


8Si/2 


11858 


1.96 [6] 


2.06 


-2] 


3.22 


0] 


54/2 


5/7/2 


11342 


5. 90 [6] 


1.52 


-1] 


3.40| 


1 1 

ij 


6P3/2 


7(^3/2 


10197 


7.61[5] 


1.19 


-2] 


1.59 


0] 


7Sl/2 


7pi/2 


39370 


2.42[6] 


5.62 


-1] 


1.46[ 


ol 

2J 


6p3/2 


7rf5/2 


10190 


4. 54 [6] 


1.06 


-1] 


1.42 


1] 


7Sl/2 


7P3/2 


37750 


2.o3[6J 


1.12 


0] 


2.79[ 


Ol 

2J 


6p3/2 


8(^3/2 


8889 


4.31 [5] 


5.11 


-3] 


5.98 


-1] 


7Sl/2 


8pi/2 


18238 


3.95[5] 


1.97 


-2] 


2.37| 


nl 

OJ 


6533/2 


8c?5/2 


8886 


2. 58 [6] 


4.58 


-2] 


5.36 


0] 


7Sl/2 


8p3/2 


18060 


5.16[5] 


5.04 


-2] 


6.00| 


nl 

OJ 


7pi/2 


6^3/2 


229885 


2.00[4J 


3.17 


-1] 


4.80 


2] 


64/2 


4/5/2 


Ibobbi 


1.11[4] 


1.77 


-1] 


6.23[ 


ol 

2J 


7Pl/2 


8Sl/2 


56883 


8. 36 [5] 


4.06 


-1] 


1.52 


2] 


64/2 


8pi/2 


39872 


3. 44 [5] 


4.10 


-2] 


2.15| 


1 1 

IJ 


7}3l/2 


7d3/2 


31928 


1.33 [6] 


4.05 


-1] 


8.51 


1] 


64/2 


8p3/2 


39032 


2.82[4J 


C A A 

6.44 


-3] 


3.31| 


nl 

OJ 


7pi/2 


8rf3/2 


21858 


8. 39 [5] 


1.20 


-1] 


1.73 


1] 


64/2 


5/5/2 


34662 


1.94 [6] 


5.24 


-1] 


2.39[ 


ol 

2J 


7p3/2 


60(3/2 


306748 


1.68 [3] 


2.37 


-2] 


9.59 


1] 


64/2 


4/5/2 


274725 


7.24[2J 


8.19 


-3] 


4.45[ 


1 1 

IJ 


7p3/2 


6(^5/2 


290736 


1 1 Of/ll 

1.12[4J 


2.22 


-1] 


8.67 


2] 


64/2 


4/7/2 


274725 


1.09 [4] 


1.64 


-1] 


8.89[ 


Ol 

2J 


7p3/2 


8Sl/2 


60643 


1.53 [6] 


4.20 


-1] 


3.36 


2] 


64/2 


8p3/2 


39200 


2.61 [5] 


4.01 


-2] 


3.10[ 


IJ 


7p3/2 


7d3/2 


33080 


2.87[5] 


4.72 


-2] 


2.05| 


1] 


64/2 


5/5/2 


34795 


1.40[5] 


2.54 


-2] 


1.74[ 


11 

J 


7p3/2 


74/2 


33003 


1.70 [6] 


4.16 


-1] 


1.81| 


2] 


64/2 


5/7/2 


34795 


2.10[6] 


5.07 


-1] 


3.49[ 


2J 


7p3/2 


8(^3/2 


22391 


1.75[5] 


1.31 


-2] 


3.87| 


0] 


8Si/2 


8pi/2 


84388 


6.59[5] 


7.04 


-1] 


3.91| 


2J 


7p3/2 


8(^5/2 


22371 


1.04[6] 


1.17 


-1] 


3.44| 


1] 


8Si/2 


8p3/2 


80710 


7.18[5] 


1.40 


0] 


7.45[ 


2J 


4/6/2 


7d3/2 


43066 


1.77[5] 


3.28 


-2] 


2.79| 


1] 


74/2 


5/5/2 


531915 


4.89[3] 


3.11 


-1] 


2.18[ 


3J 


4/6/2 


74/2 


42937 


8.24[3] 


2.28 


-3] 


1.93| 


0] 


74/2 


5/5/2 


552486 


3.12[2] 


1.43 


-2] 


1.56[ 


2J 


4/6/2 


84/2 


26560 


7.12[4] 


5.02 


-3] 


2.63| 


0] 


74/2 


5/7/2 


552486 


4.68[3] 


2.85 


-1] 


3.11| 


3J 
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noted above, the second-order E^'^^ is included in the E^^ 
value. Therefore, we use our high-precision calculation 
of E^'^'> described above to account for the contributions 
of the higher partial waves by replacement E^^^^ [^max = 
6] value with the final high-precision second-order value 

p(2) . 

£;|°, = i?s° + 41-i?(^)[W = 6]. 

The size of this correction varies from ^ 200 cm^^ for the 
lowest valence states to ~ 1 - 20 cm^^ for other valence 
states considered in this work. 

A lower number of partial waves, Zmax =6, is used also 
in the third-order calculation. Since the asymptotic l- 
dependence of the second- and third-order energies are 
similar (both fall off as we use the second-order 

remainder to estimate the numerical uncertainties in the 
third-order and in all-order corrections. 

In our calculations of the Breit contribution, we use 
the whole operator © in the first-order correction B^^\ 
while the second-order Coulomb-Breit energies are 
evaluated using the unretarded Breit operator, also 
known as Gaunt (it is described by the first term in (O). 
Usually Gaunt part strongly dominates in the Breit cor- 
rections to the valence energies ^] . Table |T] shows that 
there is strong cancelation between first and second or- 
der corrections. It is in agreement with the well known 
observation that Breit interaction for valence electrons is 
screened by the core [4l|, \^ . 

We have also estimated Lamb shift correction to va- 
lence energies. The vacuum-polarization was calculated 
in the Uehling approximation. The self-energy contri- 
bution is estimated for the s, Pi/2 and orbitals by 
interpolating the values obtained by Mohr [H, IH, |45| 
using Coulomb wave functions. We found, as expected, 
that Lamb shift correction is very small (-E ls ^ 3 cm~^ 
for In I and E ls < W cm'^ for Sn II ). This is weU 
below the accuracy of the present theory, and we neglect 
this contribution in Table HI 

Comparison of the differences 6E^'^^ = E^^^^ — i?NiST 
and SE'-^'> = E''^^^ — E'nist given in Table |T] shows that 
convergence of MBPT series is not very good for both 
In and Sn+. In particular, the second-order results for 
d-wave in In and /-wave in Sn"*" are even better than 
the third-order ones. All-order results are more accurate 
than the third-order ones, but the difference is not very 
large. For p-waves, SD calculation without the third- 
order correction overestimates valence binding energies 
and underestimates it when this correction is included. 
For d-wave, both variants lead to underestimation of the 

(3) 

binding energy and term ii'oxtra worsens the agreement 
with the experiment. 

We conclude that all-order calculation is generally 
more accurate than the third-order MBPT calculation. 
Account of the missing third-order terms does not lead 
to improvement of the accuracy. On the other hand, this 
term is generally on the order of our final difference with 
experiment and can serve as an estimate of the latter. For 



most levels, our final accuracy is better than 1%, but the 
accuracy for the d-wave of In is noticeably worse. That 
can be explained by the existence of the low-lying con- 
figuration 5s5p^ which strongly interacts with configura- 
tions 5s^nd. To account for this interaction effectively, 
one needs to consider In as a three electron atom |4(|]. 
The same reason explains mentioned above divergence 
of the SD equations for d-wave of Sn II. Interaction be- 
tween configurations 5s5p^ and 5s^ns is weaker and SD 
equations for s-wave converge for both atoms considered 
here. In the opposite parity there is no such a low- lying 
excitation of the shell 5s, so MBPT works better and no 
problems with convergence occur. 

In order to study the relative role of the valence correla- 
tions we have performed the second-order RMBPT calcu- 
lations of atomic properties of In I and Sn II considering 
these atoms as three-electron systems. Corresponding 
variant of RMBPT was developed in [13, IH, |4i,M IH ■ 
The energies of the [Re]2s^2p, [He]2s2p^ and [He]2p3 
states of B-like systems were presented in Ref. [43] . The 
second-order RMBPT was used by Johnson et al. 
to calculate [Ne]3s^3Z and [Ne]3p^3s states in Al I and 
[Xe]4/"5di°6s26pZ and [Xe]4/"5di06s6p2 states in Tl I. 
Comparing results obtained for neutral B I, Al I, and 
Tl I, we find that the discrepancy between RMBPT and 
experimental results increases significantly from B I to 
Tl I. For example, the RMBPT and NIST values of the 
ns^np ^Pj,j2 ~ '^Pi/2] splitting in cm^^ for n = 2 are 
equal to 17 and 15; for n — Z corresponding values are 
123 and 112; finaly, for n = 6 we get 6710 and 7793 
respectively. It is evident that for a light system, such 
as B I, the second-order three-electron RMBPT treat- 
ment works much better than for a heavy system, such 
as Tl I. For the latter case it is more appropriate to con- 
sider Tl I as one-electron system with [Xe]4/^'*5d^*'6s^ 
core but treat correlation more completely. It was found 
in [2l| that in such approach the discrepancy between 
the SD and NIST values of the Gs^Gp [^^3/2- '^Pi/2] split- 
ting is only 41 cm^^ instead of 1083 cm^^ obtained in 
[sa f. Alternatively, one can use CI-I-MBPT method [23j . 
where the discrepancy is 43 cm~^ [i^ . 

The main difference between configurations 
[Ni]4s24p'54di05s2nZ of In-like ions and [Ne]3s^nZ 
configurations of Al-like ions is the necessary size of the 
model space for valence electrons. For hi electrons in 
In-like ions we could not construct sufficiently complete 
three-electron model space as we did for Zl electrons. 
Additionally, in In-like ions the n — A core shell is not 
filled. Obviously, we can not exp ect the same accuracy 
as in the case of Al-like ions jSOl Is^ l . 

We tried two model spaces to evaluate energies of In- 
like ions. Firstly we constructed the model space includ- 
ing 5s, bp, and 5d electrons, [spd] model space. Secondly, 
the odd-parity model space was [5s^5p 4- 5^"^] and even- 
parity model space was [5s^5d -I- 5s5p^]. We found that 
in the second case the RMBPT energies were in better 
agreement with NIST data ^| than in the case of more 
complete [spd] model space. Theoretical values of the 
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5s^5p P-P3/2 — '^Pi/2\ splitting were equal to 2669 cm~^ 
and 4889 cm~^ in In I and Sn II, respectively. Com- 
parison of these values with the values from Ta- 
ble [J (2158 cm-i in In I and 4198 cm-^ in Sn II) shows 
that the one-electron representation with all-order treat- 
ment of correlation correction gives the results that are in 
substantially better agreement with experiment than the 
three-electron model space theory. Because of that, we 
decided not include three-electron results in the present 
paper. 



III. ELECTRIC-DIPOLE MATRIX ELEMENTS, 
OSCILLATOR STRENGTHS, TRANSITION 
RATES, AND LIFETIMES IN In I AND Sn II 

The one-body matrix element of the operator Z is 
given by [Slj : 



(10) 



where '^v,w are exact wave functions for the many-body 
"no-pair" Hamiltonian H 



H\^,) = E\^,) 



(11) 



In MBPT, we expand the many-electron function 4'^, in 
powers of Vi as 



(2) 



(3) 



(12) 



TABLE III: Wavelengths A (A), transition rates Ar (cm~^), 
oscillator strengths (/), and hne strengths 5* (a.u.) for tran- 
sitions in Sn II calculated using all-order method. 



Transition 


A 




f 
J 




s 




5pi/2 


6Sl/2 


1780 


3.17[8] 


1.47 


-11 


1.70 


nl 


5pi/2 


7Sl/2 


1170 


8.18[7] 


1.65 


-21 


1.26 


-11 


5pi/2 


8Sl/2 


1024 


7.29[6] 


1.13 


-31 


7.55 


-31 


5p3/2 


6Sl/2 


1924 


5. 76 [8] 


1.56 


-11 

-■-1 


3.90 


nl 


5p3/2 


7Sl/2 


1231 


1.35 [8] 


1.50 


-21 


2.41 


-11 


5p3/2 


8Sl/2 


1070 


1 73f7l 


1.46 


-31 


2.04 


-21 


6pi/2 


7Sl/2 


6859 


3.87[7] 


2.65 


-11 

-■-1 


1.18 


11 


6pi/2 


8Sl/2 


3735 




3.00 


-21 


7.34 


-11 


6p3/2 


7Sl/2 


7300 


7.38[7] 


2.86 


-11 

-■-1 


2.71 


11 


6p3/2 


8Sl/2 


3862 


2. 66 [7] 


2.94 


-21 


1.49 


nl 


7pi/2 


8Sl/2 


15654 


-L.UI [(J 


3.79 


-11 


3.84 


11 


7p3/2 


8Si/2 


16592 


2. 04 [7] 


4.06 


-11 

-■-1 


8.72 


11 


6Sl/2 


6pi/2 


6813 


5. 89 [7] 


4.14 


-11 


1.86 


11 


6Sl/2 


6P3/2 


6428 


6.91[7] 


8.63 


-11 


3.67 


11 


6Si/2 


7pi/2 


2852 


1.76 [5] 


2.15 


-41 


4.05 


-31 


6Si/2 


7j33/2 


2823 


1.30 [6] 


3.11 


-31 


5.78 


-21 


6Si/2 


8pi/2 


2255 


2.52[5] 


1.93 


-4] 


2.86 


-3] 


6Si/2 


8p3/2 


2245 


1.38[1] 


2.08 


-8] 


3.08 


-7] 


7Si/2 


7pi/2 


17218 


1.21[7] 


5.71 


-1] 


6.69 


1] 


7Si/2 


7p3/2 


16210 


1.41[7] 


1.18 


0] 


1.30 


2] 


7Si/2 


8pi/2 


6625 


4.04[5] 


2.72 


-3] 


1.20 


-1] 


7Si/2 


8p3/2 


6545 


9.67[5] 


1.27 


-2] 


5.54 


-1] 


8Si/2 


8pi/2 


34495 


3. 62 [6] 


6.95 


-1] 


1.64 


2] 


8Si/2 


8p3/2 


32425 


4.26[6] 


1.43 


0] 


3.16 


2] 



The denominator in Eq. ()10|) arises from the normaliza- 
tion condition that starts to contribute in the third order 
[53I. In the lowest order, we find 



(13) 



where Zwv is the corresponding one-electron matrix ele- 
ment. Since '^w^ is a DF function we use Z^^ designation 
instead of Z^^^ below. 

The second-order Coulomb correction to the transition 
matrix element in the case of V^^^ DF potential is given 

by m 



E 



gn 



wavn Qwanv )'^na 



Sn, ~r £"71 



(14) 

The second-order Breit corrections are obtained from 
Eq. ([T4|) by changing gijki to bijki [1^. The third-order 
Coulomb correction is obtained from Eqs. (|10p and ([l3 
as 

= (vI/(„2)|^|^(0)^^^^(0)|^|^(2)^_^^^a)|^|vI,W 



7(1) r 



(vI/W|vI/W}-|-(vI;a)|vi,a)) 



(15) 



where the last term arises from the normalization condi- 
tion. In Ref. [5§|, contributions to Zwl were presented 
in a following form: 



NORM 



(16) 



The first term here corresponds to the well known ran- 
dom phase approximation (RPA). Though RPA corre- 
sponds to the summation of certain MBPT terms to all 
orders, it is possible to include it here using the procedure 
described in Ref. [Ill . Next term Z^*-* corresponds to the 
correction which arise from substituting DF orbitals with 
Brueckner ones. The last two terms in Eq. describe 
structural radiation, Z^^, and normalization, Z^"^^, 
corrections. 

In the all-order SD calculation, we substitute the all- 
order SD wave function into the matrix element ex- 
pression given by Eq. ([TU)) [3l|: 



wv 



^{l + N^){l + N,) 



(17) 



i,^^ is the DF matrix element p3p and the terms 
k — a- ■ -t are linear or quadratic function of the ex- 
citation coefficients introduced in Eq. ([5]). Normalization 



where 2,, 
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terms Ny^jj^ are quadratic functions of the excitation co- 
efficients. This expression completely incorporates Z^'^^ 
and certain sets of MBPT terms are summed to all or- 
ders 'si'l . The part of the fourth-order correction that is 
not included in the SD matrix element p7D was recently 
discussed by Derevianko and Emmons [55], but we do 
not include it here. 

In Tables [TT] and IIIIl we present theoretical transition 
rates Ar, oscillator strengths /, and line strengths S for 
El transitions between low-lying states of In I and Sn II, 
respectively. These results are obtained by combining 
all-order El amplitudes pT|) in the length gauge and 
theoretical energies Ef^^ from Table U using well-known 
expressions (see, for example, Ref. 25]). 

Calculation of the transition amplitudes provides an- 
other test of the quality of atomic-structure calculations 
and another measure of the size of the correlation correc- 
tions. In Tables [TVl and fVl we compare our results with 
available experimental data. For convenience, we also 
present theoretical and experimental wavelengths for all 
transitions. There is good agreement with experimental 
results for the strongest lines of In. For Sn II, agreement 
is also good with exception of the 6s - 6^3/2 transition 
where experimental value is much smaller than the cal- 
culated one. Note that the theory and experiment are in 
good agreement for the 6s - 6pi/2 transition. 

We also use El transition rates to calculate the life- 
times of low- lying levels of In I and Sn II. We compare 
these lifetimes r^^°' with available experimental mea- 
surements in Table IVIl For Tdj-levels, the measure- 
ments from Refs. 8, ^, JJJ gave rather different lifetimes. 
Our calculations support the shorter times obtained in 
Ref. a. 



IV. HYPERFINE CONSTANTS FOR INDIUM 

Calculations of hyperfine constants follow the same 
pattern as calculations of El amplitudes, described in 
the previous section. The value of the nuclear magnetic 
moment for "'^^^In used here is taken from fsTj. Hyperfine 
constants for another odd isotope, ^^"^In, can be obtained 
using the scaling factor 0.99785, which is indistinguish- 
able from unity within the accuracy of the present theory. 
In contrast with dipole amplitudes considered above, the 
hyperfine structure is sensitive to the wave function at 
short distances and to very different types of correlation 
corrections. 

Table lVlIl shows that SD method significantly improves 
DF values of the hyperfine constants of the lowermost lev- 
els. It is rather unusual that correlation correction to the 
hyperfine structure constant of 5p3/2 level is so small. For 
other p3/2-levels, correlation corrections are comparable 
to the initial DF contribution. This situation is more typ- 
ical for other atoms with ns^np3/2 configuration, such as 

Ti gi,!!!. 



TABLE IV: Oscillator strengths / and wavelengths A (A) in 
In I. The SD data (f^^) are compared with semi-empirical 
calculationsX/^^) from Ref. [it] and experimental data (f^^^) 
from Ref. 

Lower Upper A^'p'- f° /^'p'' 



V. CONCLUSION 

Summarizing results of the previous sections we can 
make several conclusions. We have seen that all-order SD 
calculations, when converge, provide an improvement to 
the third order MBPT calculation. Convergence of the 
SD-equations is hampered by the existence of low-lying 
excitations from the uppermost core shell 5s. The lowest 



5?7i /o 


6Sl/2 


4153 


4102 


^V-i 10 


6Si/2 


4576 


4511 


5pi/2 


7tj , 
' Sl/2 


97Q9 


A i 04 






2977 


2933 


5pi/2 


8Si/2 


2493 


2460 


5p3/2 


8S1/2 


zDoy 


ZDUZ 


5pi/2 


9Si/2 


2370 


2340 


5p3/2 


9Si/2 


2502 


2468 


5pi/2 


5^3/2 


3045 


3039 


5p3/2 


5^3/2 


3266 


3259 


5p3/2 


5^5/2 


3264 


3256 


5pi/2 


6^3/2 


2578 


2560 


5p3/2 


6^3/2 


2734 


2713 


5p3/2 


6^5/2 


2734 


2710 


5pi/2 


7^3/2 


2410 


2388 


5p3/2 


7^3/2 


2547 


2523 


5p3/2 


7^5/2 


2546 


2521 


5pi/2 


8^3/2 


2329 


2306 


5p3/2 


8^3/2 


2456 


2432 


5p3/2 


805/2 


Z40D 


z4oy 


6z?i /o 


7Sl/2 


22594 




6p3/2 


7Sl/2 


24184 






8Sl/2 


11463 






8Si/2 


11858 




6pi/2 


9Si/2 


9264 




6p3/2 


9Si/2 


9520 




6Sl/2 


6pi/2 


13669 




6Si/2 


6P3/2 


13146 




6Si/2 


7pi/2 


7002 




6Si/2 


7p3/2 


6949 




6Si/2 


8pi/2 


5806 




6Si/2 


8p3/2 


5787 





0.133 


0.137 


0.14 


0.142 


0.153 


0.15 


0.016 


0.0158 


0.017 


0.015 


0.161 


0.017 


0.0054 


0.00541 


0.006 


0.0051 


n r\c\ on 

0.00539 


0.006 


0.0025 


0.00256 


0.0029 


0.0024 


0.00254 


0.0026 


0.361 


0.51 


0.36 


0.040 


0.056 


0.06 


0.353 


0.49 


0.37 


0.072 


0.11 


0.043 


0.0077 


0.011 


0.006 


0.068 


0.10 


0.052 


0.026 


0.039 


0.006 


0.0027 


0.0033 


0.0014 


0.024 


0.035 


0.009 


0.012 


0.017 


0.0003 


0.0013 


0.0016 




A1 1 

U.Oii 


O.uio 


U.UUio 


0.266 


0.274 




0.279 


0.287 




0.0222 


0.233 




U.UzU ( 


n 01 Q 
U.zlo 




0.00729 


0.00764 




0.00664 


0.00702 




0.402 


0.467 




0.813 


0.944 




0.0103 


0.0110 




0.0284 


0.0207 




0.00206 


0.00223 




0.00642 


0.00704 
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TABLE V: Transition probabilities A (in lO'' s~^) and wave- 
lengths A (A) in Sn II. Our SD results are compared with 
experimental data from Ref. f^. 



Lower 


Upper 


A^D 


^cxpt 


^SD 




6Sl/2 


6Pl/2 


6813 


6844 


5.89 


5.8±1.1 


6Pl/2 


7Sl/2 


6859 


6761 


3.87 


4.2±0.1 


6Si/2 


6P3/2 


6428 


6453 


6.91 


5.2±1.0 


6P3/2 


8Si/2 


3862 


3841 


2.66 


2.5±0.5 


6Pl/2 


8Sl/2 


3735 


3715 


1.45 


1.8±0.3 



TABLE VI: Lifetimes r in ns for the nl levels in indium. The 
SD data are compared with experimental results from (a)- 
Ref. H, (6)-Ref. ^, (c)-Ref. fTll|, and (d)-Ref. 



Level 






Level 




6Si/2 


7.04 


7.5±0.7" 






7Sl/2 


21.5 


19.5±1.5'=; 19.5±1.5''; 27±6'' 


6pi/2 


69.7 


8Sl/2 


47.7 


53±5'=; 55±6'' 


7Pl/2 


219. 


9Sl/2 


89.4 


118±10'=; 1Q4±12'' 


7p3/2 


192. 


5^3/2 


6.45 


6.3±Q.5'' 


8pi/2 


473. 


5^5/2 


6.78 


7.6±0.5'' 


8p3/2 


414. 


6^3/2 


19.2 


21±3" 


4/5/2 


70.4 


6^5/2 


20.1 


22±3"; 18.6±1.5'=; 18.6±1.5'' 


4,f7/2 


70.4 


7^3/2 


42.0 


50^:5"; 200±4'' 


5/5/2 


125. 


7d5/2 


44.0 


50±5''; 154±10'=; 147±10'' 


5/7/2 


125. 


8d3/2 


75.7 


317±22'= 






8d5/2 


77.2 


300±60'=; 238±20'' 






6P3/2 


63.7 


55.0±4'* 







such excitation corresponds to configuration 5s5p^ that 
has positive parity and primarily affects SD-equations for 
the valence d-wave. Because of that, we were not able to 
solve these equations for Sn II. To avoid this problem 
one has to exclude 5s-electrons from the core and con- 
sider In I and Sn II as three-electron systems. However, 
our attempt to treat In I and Sn II as three-electron sys- 
tems within second-order RMBPT for the valence model 
space, as suggested in Refs. [i^, resulted in rather poor 
agreement with experimental spectra. We conclude that 
valence correlations for atoms in question can not be 
accurately accounted within model space approach. It 
would be interesting to perform CI+MBPT calculations 
[2^, but it goes beyond the scope of the present paper. 



Another interesting observation concerns the addition 

(3) 

of the missing part of the third-order term -Eextra to the 
SD results. It was suggested in Ref. [s^l to add this 
term, so that all third-order terms are accounted for. For 
heavy alkali-metal atoms omission of this term leads to 
significant discrepancies of the all-order values with ex- 
periment. One can see from Table HI that for atoms con- 
sidered here this term does not improve agreement with 
experimental energies. We have also found that first- and 

TABLE VII: Hyperfine constants, A (in MHz) for "^In 
(1=9/2, ^i=5.5408 [131). Dirac-Fock (DF) and all-order 
(SD) calculations are compared to experimental values from 
Ref. Il5i - (a), Ref. [H - (6), and Ref. H - (c). 



Level 


DF 


SD 


Exper. 


Level 


DF 


SD 


6si/2 


983.0 


1812 


1685" 


5d3/2 


4.365 


-11.48 


7Sl/2 


335.6 


544.5 




6^3/2 


2.305 


-11.20 


8si/2 


153.6 


240.8 




7^3/2 


1.275 


-7.692 


9Sl/2 


83.10 


128.1 




8(^3/2 


0.805 


-5.385 


5pi/2 


1780 


2306 


2282** 


5^5/2 


1.862 


47.83 


6pi/2 


222.7 


263.2 




60(5/2 


0.981 


30.81 


7pi/2 


85.15 


95.61 




7^5/2 


0.543 


18.95 


8pi/2 


41.90 


45.97 




8d5/2 


0.342 


12.59 


5P3/2 


267.8 


262.4 


242.2'= 


4/5/2 


0.0611 


0.1871 


6P3/2 


35.69 


77.82 




5/5/2 


0.0316 


0.1055 


7P3/2 


13.71 


30.83 




4/7/2 


0.0339 


0.2293 


8p3/2 


6.767 


15.42 




5/7/2 


0.0176 


0.1658 



second-order Breit corrections tend to cancel each other 
in agreement with [i^; final Breit corrections are small 
and can be neglected within present accuracy of the the- 
ory. 
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